New precision measurement of the J/ip- and ^'-meson masses* 
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A new high precision measurement of the J/ip- and f//-meson masses has been performed at 
the VEPP-4M collider using the KEDR detector. The resonant depolarization method has been 
employed for the absolute calibration of the beam energy. The following mass values have been 
obtained: 

M J/V , = 3096.917 ± 0.010 ± 0.007 MeV, 

M^/ = 3686.111 ±0.025 ±0.009 MeV. 

The relative measurement accuracy has reached 4 • 10~ 6 for J/tp and 7 • 10~ 6 for ip', approximately 
3 times better than in the previous precise experiments. 
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I. INTRODUCTION 

This work continues a series of experiments on the pre- 
cise determination of the onium resonance masses at the 
electron-positron collider VEPP-4: J/tp , tp' (OLYA de- 
tector) and T, T', T" (MD-1 detector) 0,111111. 
A few years ago the values of the masses obtained have 
been rescaled |(| to take into account the progress in 
the electron mass measurements H |tj . 

VEPP-4 experiments employed the resonant depolar- 
ization method [Tol[Tl| for the absolute beam energy cali- 
bration and achieved the relative mass accuracy of 1- 10 -5 
for the T-family and of 3T0 -5 for the ^-family. The reso- 
nant depolarization experiments on bottomonium masses 
were also performed with the CUSB detector at CESR 
P3 (T) and with the ARGUS detector at DORIS 
(T'). The accuracy of the J/-0-mass measurement was 
improved in the Fermilab pp-experiment E760 |14| to 
1.2 ■ 1CU 5 using the ip' mass value from Ref. pj. 

The goals of this work were to further improve the ac- 
curacy of the J/ip- and -^'-masses and develop the reso- 
nant depolarization technique at the upgraded VEPP-4M 
collider for future experiments. 

The first precise measurement of the J /ip and ip' me- 
son masses [l| set the mass scale in the range around 
3 GeV which provided a basis for the accurate determi- 
nation of the charmonium state location. At present the 
charm meson family is a good test bench for QCD and 
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quark potential models predictions in which masses of 
the open and hidden charm can be calculated with good 
accuracy. Another fundamental application of the men- 
tioned measurements is the r-lepton mass determination 

M 

Substantial improvement in the beam energy accuracy 
obtained by the presented experiment sets a new stan- 
dard of the mass scale in the charmonium range. 



II. BEAM ENERGY DETERMINATION 
TECHNIQUE 

A. Resonant depolarization method 

Electrons and positrons in storage rings can become 
polarized due to emission of synchrotron radiation ac- 
cording to the Sokolov-Ternov effect 0] . Spins of polar- 
ized electrons precess around the vertical guiding mag- 
netic field with the precession frequency f2, which in the 
plane orbit approximation is directly related to the par- 
ticle energy E and the beam revolution frequency lo: 



n/uj = i ± 7 • ///mo = i + v 



(1) 



where 7 =E/m , m c is the electron mass, // and no are 
the anomalous and normal parts of the electron magnetic 
moment. The v is a spin tune, which represents the spin 
precession frequency in the coordinate basis related to 
the particle velocity vector. 

The precession frequency can be determined using the 
resonant depolarization. To this end one needs a polar- 
ized beam in the storage ring which is affected by the ex- 
ternal electromagnetic field with the frequency f2 d given 
by the relation 



(2) 
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with any integer n (for VEPP-4M in the J/ip region 
n = 3). 

The precession frequency is measured at the moment of 
the polarization destruction detected by the polarimeter, 
while the depolarizer frequency is being scanned. The 
process of forced depolarization is slow enough compared 
to the period of the synchrotron oscillations of the par- 
ticle energy. This allows to determine the average spin 
tune {u) and corresponding average energy of the par- 
ticles (E) with higher accuracy than the beam energy 
spread o~e- 

Due to modulation of the precession frequency by par- 
ticle orbital motion, the resonant depolarization could 
happen at the sideband resonances, which are distant 
from the main one by multiples of the synchrotron and 
betatron frequencies. Besides, it could happen at the 
weak sideband resonances caused by extraneous low fre- 
quency modulation of the guide field, caused for example 
by pulsations in the power supply system (50 Hz, the en- 
ergy shift of about 25 keV). Therefore, it is necessary to 
identify the main resonance by special means. 

It should be noted that the average energy of the 
beam particles (E) differs from the energy of the equi- 
librium particle E s because of the radial betatron oscil- 
lations. The effect is proportional to the betatron ampli- 
tude squared and is mainly due to the nonlinearity of the 
guide field. It also determines the spin resonance natu- 
ral bandwidth 0|- I n this experiment the observed full 
bandwidth was about 5 keV in beam energy units. 

Formula (JJJ gives the value of 7 averaged over the 
beam revolution time. Thus, for a symmetric machine, 
it corresponds to the energy in the interaction point. 

The method described has been developed in Novosi- 
birsk and first applied to the </>-meson mass measurement 
at the VEPP-2M storage ring ^(J- The comprehensive 
review of the resonant depolarization technique and its 
applications for particle mass measurements can be found 
in [13. 



B. VEPP-4M polarimeter 



4 Hz band. This provides controllable conditions [2lJ for 
the depolarization at the main spin resonance at the min- 
imal level of the depolarizer power, which corresponds to 
~ 2 • 10~ 6 rad spin rotation per a single pass of the par- 
ticle, and with the frequency-tuning rate of 0.2 Hz/sec. 

The characteristic jump in the relative rate of scat- 
tered electrons at the moment of resonant depolarization 
is 3.0 ~ 3.5% with the statistical error of 0.3 ~ 0.4% for 
the beam polarization degree higher than 50%. Typical 
behavior of the rate ratio is shown in Fig. The linear 
growth before the depolarization reflects the difference 
in the bunch life times due to polarization dependence of 
the intra-beam scattering cross section. 
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FIG. 1: The variation of the coincidence rate ratio for polar- 
ized and unpolarized beams during the energy calibration (Ei 
is the energy calculated from the magnet currents, Erede 
is the energy determined by resonant depolarization, Fd ep is 
the corresponding depolarizer frequency; the vertical numbers 
show the instant depolarizer frequency). 

The characteristic uncertainty of the beam energy cal- 
ibration due to the depolarization procedure is 1.5 keV 
(for more detail see Ref. [2fjj1- 



The polarimeter unit is installed in the technical 
straight section of VEPP-4M and consists of the po- 
larimeter employing the spin dependence of the intra- 
beam scattering (Touschek) effect [l^ and TEM wave- 
based depolarizer [20| . 

The polarimeter detects Touschek electron pairs with 
the help of two movable scintillation counters placed in- 
side the beam pipe pockets. We use the "two bunches" 
compensation technique, in which relative rates of scat- 
tered particles from unpolarized and polarized beams are 
compared. 

The rate of Touschek electrons is 3 -j- 12 kHz at the 
beam current of 2 4 mA. The depolarizer frequency is 
scanned with a step of 2 Hz by the computer controlled 
synthesizer with the intrinsic bandwidth of ~ 1 Hz. How- 
ever, the frequency line is artificially broadened up to the 



C. Accuracy of single energy calibration 

The achievable accuracy of the J/ 1\) and ip' mass mea- 
surement was analyzed in |22| . Since that time, under- 
standing of some systematic effects improved. In partic- 
ular, the energy shift due to the vertical closed orbit dis- 
turbances turned out to be much less than was expected 
in |13 and |2^ . 

The relation is broken in the radial magnetic field 
and the vertical electric field (used for the electrostatic 
beam separation) because of the orbit non-planarity. The 
influence of these fields in the second order of perturba- 
tion theory can be expressed in terms of Fourier har- 
monics of the vertical closed orbit disturbances. The 
accurate analysis of the effect and the numerical simu- 
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lation gave the energy shift of — 0.4 ± 0.3 keV for J/ip 
and -0.3 ±0.2 keV for ip' . 

The uncertainty estimates for the mean value of the 
beam energy in the interaction point for a single calibra- 
tion are collected in Table [I] 

The energy value obtained in a single calibration is 
biased due to a non-zero spin resonance width. The re- 
quired correction (w 2.5 keV) can be determined with 
accuracy better than 1 keV using a few calibrations with 
opposite directions of the depolarizer frequency scan. 

The coherent energy loss asymmetry arises from the 
difference of the resistive impedance of the beam pipe in 
two half-rings. As a result, the energy of electrons and 
positrons in the interaction point (I. P.) differs from the 
energy value obtained by the resonant depolarization. 

The uncertainties not exceeding 0.1 keV are not shown 
in the table including those due to the non-zero mo- 
mentum compaction factor and the longitudinal magnetic 
fields 22]. 

The uncertainty of a single energy measurement does 
not directly contribute to the systematic error of the me- 
son mass. Thus, the effect of the energy loss asymme- 
try in the half-rings has an opposite sign for e + and e~ 
and cancels in the linear approximation. The contribu- 
tion of the precession and revolution frequency measure- 
ments has mainly statistical nature and becomes negli- 
gible when a large number of calibrations is used. At 
the same conditions the uncertainty due to the non-zero 
spin resonance width vanishes provided that frequency 
scan directions alternate. 

On the other hand, new sources of the systematic er- 
ror come into play when a long-term experiment with 
colliding beams is considered. This is discussed below in 
Sec. lIVBl Sec. Ivland. finally, Sec. IVIlil when the essen- 
tial features of the experiment are described. 
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FIG. 2: The data acquisition scenario for J/ip (circles) and 
the actual points of the second scan (points with error bars). 
The solid line shows the second scan fit. 

and 7 the required integrated luminosity was reduced by 
the factor of 2. Such a 5 + 2 scheme does not minimize 
statistical and systematic errors for a given experiment 
duration but allows one to apply the \ 2 criterion to the 
results of a single scan. 

Before data acquisition at point 1, the beam energy 
calibration was made to fix the current energy scale. At 
points 2-6 the calibrations before and after data tak- 
ing were performed with the opposite direction of the 
depolarizer frequency scan. The point 7 requires no en- 
ergy calibration. From 2 to 5 ring fillings were neces- 
sary at each point to collect the wanted integrated lu- 
minosity. The injection from VEPP-3 occurred at the 
required set-up energy without any intentional change of 
the VEPP-4M magnet currents. 

On completion of the scan, the VEPP-4M magnetiza- 
tion cycle was performed and the whole procedure was 
repeated. 

The set-up parameters of the collider and the results of 
the current, magnetic held, temperature and orbit mea- 
surements have been stored in the database. 



III. EXPERIMENT DESCRIPTION 

The first part of the experiment consisted of three 
scans of the J/ip -region (the integrated luminosity 
J Ldt w 40 nb _1 , the beam energy spread erg » 0.6 MeV) 
and three scans of the -(//-region (J Ldt » 76 nb _1 , 
o~ e ~ 0.9 MeV). Then the betatron and synchrotron 
damping decrements of VEPP-4M were rearranged to re- 
duce the energy spread to 0.45 MeV and the fourth scan 
of J /if) was performed (/ Ldt w fO nb _1 ). The goal was 
to verify possible systematic errors related to the collider 
operating mode and the beam energy spread. 

The beam polarization time in the VEPP-4M ring is 
about 100 hours at the ^-energy region. For the energy 
calibration runs, the beam spent the time sufficient for 
the polarization in the booster ring VEPP-3 (2.5 hours 
at J/ip and about 1 hour at ip') and was injected to 
VEPP-4M without essential loss of the polarization de- 
gree. 

During the scan the data were collected at 7 points 
of the resonance excitation curve (Fig. [2J ■ At points 1 



IV. ASSIGNMENT OF ENERGY TO DATA 
ACQUISITION RUNS 

Unlike the experiments 0- C3- C3| ■ the energy calibration 
during the data acquisition was not possible, therefore 
it was necessary to assign the energy to each data ac- 
quisition run using the results of the energy calibration 
runs. 



A. Stability runs and energy prediction function 

To make the reliable energy assignment, two stability 
runs consisting of the continuous series of energy cali- 
brations without any intentional magnet current changes 
were performed, one after the third scan ofj/j' (March 
2002) at the beam energy E « 1846 MeV [23 and the 
other one after the fourth scan of J/tp (May 2002) at 
E w 1550.6 MeV. Their results are presented in Figs. l3l4l 

One can see a rather slow energy variation and the 
day-to-night oscillation with the amplitude growing from 
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TABLE I: Single energy calibration uncertainties in the vicinity of J/ip and tp' (keV). 

Nature J/ip ip' 



Source 

Vertical orbit disturbances 
Spin resonance width 

Coherent energy loss asymmetry 

Precession frequency measurement accuracy 

Revolution frequency measurement accuracy of 10 

Sum in quadrature 

* — correction uncertainty 



Systematic " ~ " 0.3* 0.2* 

Systematic, depends on the frequency 1.0* 1.0* 
scan direction 

Systematic, charge depending 0.6 1. 

Statistical 1.2 1.5 

Statistical 0.2 0.2 

1.7 2.1 




61 62 111 



(b) day of the experiment (b) day of the experiment 



FIG. 3: The results of the first stability run: (a) — parabolic 
fit of the energy deviation, (b) — the same with the fit sub- 
tracted (the error bars show the mean deviation from this 
fit). 



FIG. 4: The results of the second stability run: (a) — energy 
deviation, (b) — the average orbit position given by BPMs 
where available (the error bars show the mean deviation from 
the fit). 



4 keV in March to 45 keV in May. The energy oscilla- 
tions with the amplitude of (45 ± 5) keV observed since 
the 110th day of the experiment are correlated with time 
dependence of the average orbit position with the am- 
plitude of Ai? = (51 ± 10) /im. The estimate of the cor- 
responding energy deviation is (77 ± 15) keV under the 
assumption that the collider expands and shrinks uni- 
formly. 

The large and relatively fast energy variations do not 
allow us to use the mean energy of the two calibrations 
surrounding a data acquisition run as the beam energy 
for this run, as was supposed initially. Instead, energy 
prediction functions have been suggested, which employ 
the results of the field measurements in some magnets by 
nuclear magnetic resonance (NMR) and the temperature 
measurements and include the explicit time dependence 
as a substitute of variables, which were not monitored 
(the effect of the tunnel wall temperature on the ring 
perimeter etc.). The orbit measurements at VEPP-4M 
are not accurate and comprehensive enough to be used for 
the energy prediction (there are ten independent bending 
magnet power supplies, sixty radial correctors and only 
fifty four beam position monitors). 

Several prediction functions have been tried of the 



form: 

E p = 3? ■ Hnmr ■ (1 + x • (T ring — Tnmr)) x 

8E on -exp (-— ) -\-SE cyc i e -cxp (-hy^l) + (3) 

A(t) • cos - ip(t)j + E (Ai, t), 

where Hnmr is the field in the off-ring calibration mag- 
net with the temperature of Tmmr] T ring , T air and 
T W ater are the average values of the ring and cooling agent 
temperatures, t is current time, t on and t cyc i e denote time 
elapsed since the last switching on the collider and the 
last magnetization cycle, respectively. The x, 6E on , 
5E cyc i e , r on , T cycie and t are free parameters determined 
by the fit of all energy calibrations performed in a certain 
operation mode (J/ip scans I-j-III, ip' , J/ip scan IV). 

The first term is the basic parameterization of the mag- 
netic field integral behavior. The exponential terms in- 
troduce the magnetic field relaxation; ignoring this effect 
results in a bias in the parameters x etc. The oscil- 
lating term describes the variations of the ring perimeter 
under the assumption that the Fourier harmonics around 
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day of the experiment day of the experiment 



FIG. 5: The predicted energy with all aperiodic dependences 
removed (the error bars show the mean deviation from the 
prediction) . 



FIG. 6: Aperiodic energy dependence on time due to switch- 
ing on the ring and the magnetization cycles (the error bars 
show the mean deviation from the predicted energy). 



il = 2it/t dominate. The A and <p are not constant due 
to season variations in properties of the ground surround- 
ing the machine tunnel. The term Eo(Ai,t) takes into 
account the energy variation due to a few relatively large 
adjustments of the current Ai in some machine lattice 
elements. 

The prediction functions differ by the choice of f(T), 
A(t) and f{t). These (simplest) functions, as well as 
.E?o(Ai,f), have additional free parameters. 

The results of the best energy prediction for May 2002 
(the fourth scan of J/ip and the second stability run) 
are illustrated by Figs. |S] and It should be noted 
that the time dependence presented in Fig. is partially 
compensated by the temperature and field strength de- 
pendences. Fourteen free parameters were used to fit 
28 points shown in these figures. The value of the chi- 
squared per degree of freedom was employed to estimate 
the (mean) error of the energy prediction by requiring 
X 2 /NdoF — 1- According to the fit, the energy oscilla- 
tion period r = 1.02 ± 0.02 days, so we have fixed it at 
the value of 1. 

The appearance of strong oscillations (FigEl the 100th 
day of the experiment) and their further growth (between 
107th and 110th days) can be probably explained by the 
change of the thermomechanical properties of the ground 
surrounding the half-ring tunnels. These properties can 
change abruptly at the moment when the melting front 
reaches the tunnel (in Novosibirsk it occurs in May) . Ac- 
cording to the x 2 criterion, the sudden growth of the os- 
cillation amplitude on the 100th day is much more prob- 
able than the gradual one. 

The direction of the depolarizer frequency scan was 
taken into account in the prediction function fit. This 
allows us to determine the spin resonance half-width of 
2.4 ± 0.7 keV in beam energy units and makes the predic- 
tion function value unbiased. The accuracy of the energy 
prediction varies from 6 to 8 keV during the whole ex- 
periment (218 calibrations). 



B. Energy assignment accuracy 

The uncertainty of the energy prediction of (6 -j- 8) keV 
includes the statistical error of the single energy cali- 
bration and all the uncertainties due to the run-to-run 
difference. Among the sources of the latter, the radial 
closed orbit variations and the RF frequency drift dom- 
inate |23. The model of statistically independent orbit 
perturbations used in overestimates the effect of the 
radial closed orbit variations. 

The significant difference between the energy calibra- 
tion runs and the data acquisition runs was the beam 
separation at the additional interaction point. This sep- 
aration could cause some small sign-of-charge dependent 
shift of the e + - and e _ -beam energy. During the cali- 
brations, the separation was turned off, so in the linear 
approximation the energy of the single beam is equal to 
the mean energy of two beams with the separation on. 
During the data acquisition, the separation was turned 
on to provide normal conditions for beam-beam effects. 
This increases the length of the beam orbit in the straight 
section and decreases it in the half-rings causing the 
equal energy shift for the e + - and e~-beams. The es- 
timated energy shifts are (—1-7 ± 0.2) keV at J/tp and 
(-2.0 ±0.2) keV at ip' . 

Due to the collision effects, the amplitude of the ra- 
dial betatron oscillations in the data acquisition runs is 
bigger than that during the energy calibration, there- 
fore, the average particle energies are not the same (see 
Sec. Ill A|) . The corresponding systematic errors of the 
energy assignment do not exceed 1.5 keV for J/tp and 
1.8 keV for ip'. 

With the two exceptions mentioned above, the energy 
calibration runs and the data acquisition runs do not 
differ essentially and we assume the same statistical ac- 
curacy of the energy prediction for them. 

Estimates of the beam energy uncertainty for the data 
acquisition runs are presented in Table[n] The statistical 
component of the beam energy uncertainty contributes to 
the meson mass error with the factor oc 1/ VN, where N 
is the number of calibrations. The accurate calculation 
of this contribution as well as that due to the choice of 
the energy prediction function can be performed using 
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TABLE II: The energy uncertainty for the data acquisition runs in the vicinity of J/ip and ip' (keV). 



Source 


Nature 


J/V 


</>' 


Energy prediction 


Statistical 


7.6 


6.5 


Radial betatron oscillations 


Systematic 


<0.7 


<0.9 


Single energy calibration (Sec. Ill CI 


Systematic 


0.3 


0.2 




Charge depending 


0.6 


1.0 


Beam separation in the additional LP. 


Systematic 


0.2* 


0.2* 


Sum in quadrature 




^7.7 


^6.6 



* — correction uncertainty 



the actual resonance curve fit (see Sec. IVIlfll . 



V. DETERMINATION OF MEAN COLLISION 
ENERGY 

The production rate of the resonance of mass M at the 
given collider energy E is determined by the probability 
of e + e~-collisions with the invariant mass W ~ M. Tak- 
ing the angular spreads 9 X , 9 y and the energy spread cte 
into account, one has after averaging over the particle 
momenta 



{W) 



±(ei + e 2 y )E- 

E ((E + ) {E_)f 
2E 



a 



(4) 



AE 



The value of E is determined by the resonant depolar- 
ization; the last term is due to the difference of the co- 
herent energy loss in two half-rings rSec.lTTCf. For the 
VEPP-4M conditions at the ■(/'-energy region the correc- 
tion is about 0.2 keV, thus we assume W = E + + E- for 
each collision. 

For beams with the Gaussian energy spread in the pres- 
ence of the electrostatically induced vertical dispersion 
ip* and the beam impact parameter A y , the differential 
luminosity can be written as 



dL(E, W) 



1 



dW 



A-K<j%{W/2)a*JW/2) ^ aw 



exp ■ 



W -2E (T W %A y 



o~w 



2Ea 2 



Ac* 



(5) 



where fu is a revolution frequency, N + and iV_ are the 
bunch populations. The transverse beam sizes in the 
interaction point a*, a* effectively depend on the sum 
E + + E- due to the collider (3- function chromaticity (the 
formula is valid to the first order of this effect and as- 
sumes the beam symmetry). According to (JSJ), the mean 
collision energy (W) L = JwdL(E,w) ^ 2E, which leads 
to the systematic error in the resonance mass. 

The vertical dispersion \ip*\ ss 800 [im with opposite 
signs for e + and e~ appears in VEPP-4M due to the 
beam separation in the additional interaction point. The 
residual orbit perturbations related to this separation re- 
sult in the beam misalignment in the experimental LP. 



characterized by A y . It leads simultaneously to the lumi- 
nosity loss and the collision energy shift. To avoid that, 
the voltage on the experimental LP. separator plates was 
tuned to have maximum luminosity in each run with the 
accuracy better than 3% so that A y < 2.5 /im at the 
beam size a* ~ 7 /im. Thus, the uncertainty of the col- 
lision energy W in the run is less than 10 keV for J/ip 
and less than 18 keV for ip'. In the resonance mass error, 
this uncertainty is suppressed oc where N > 100 

is the number of runs (see Sec. IVIII^) . 

The /3-function chromaticity also leads to the shift of 
the mean value of the collision energy. For the given emit- 
tance e y , the vertical beam size a* — yjt y l3*. Using ap- 
proximations /3 y (W/2) « (3* (l+d\nP*/dE (W/2-E)) 
and 



Jl + din f3*/dE (W/2-E) 



exp <; --dlnf3*/dE (W/2-E) 



(6) 



in ©, one obtains 8 (W) = -1/4 d\n(3*/dE a^. The 
effect of the radial /3-function chromaticity is suppressed 
by the factor of (a* p/cr*) 2 ~ 0.1 where a* ^ is the beta- 
tron radial size and a* is the total size including the dis- 
persion. The (3*(E) measurements at VEPP-4M gave the 
shifts of -4 ± 2 keV for J/ip scans I-e-III, -1.5 ± 0.7 keV 
for J/ip scan IV and +5 ± 2.5 keV for ip' (the chro- 
maticity was partially compensated after the third scan 
oiJ/tp). 

At the 1 keV level of accuracy the potential energy 
of colliding particles inside the beams should be taken 
into account. The effective energy of the electron is 
Ekinetic + U/2, where the potential energy U is due to 
its Coulomb interaction with all other electrons of the 
beam. For the flat beam with the logarithmic accuracy 



U 



e 2 N , D 
—= — In — 



(7) 



where N is the bunch population, a z is the longitudinal 
bunch size and D is the beam pipe diameter (in the beam 
rest frame the interaction of particles at longer distances 
is screened out). The kinetic and potential energies in the 
LP. differ from those in the ring because of the difference 
in the beam and beam pipe sizes, but the total energy 
conserves during the revolution, therefore 

Ekinetic. I .P. "t - Uj_p_/2 — Ekinetic^ing "t - U r i n g /2 (8) 
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At the moment of the annihilation the total energy of the 
e + e~ pair transforms to the product mass, thus 

W =2 ■ {EuneticJ.P. + Ui.p) = 

inetic,ring + Uring + Ul.P. 

The resonant depolarization result is ss Ekinetic,ring, 
therefore the collision energy shift SW = Ui.p. + U r i ng - 
The energy losses ignored in this consideration does not 
change the final result. 

For the the actual values of the beam currents and sizes 
it leads to a correction of (2 ± 1) keV for J /ip and ip' . 

VI. EVENT SELECTION AND LUMINOSITY 
MEASUREMENTS 

A. Detector and trigger 

The KEDR detector [24[ consists of the vertex detec- 
tor, the drift chamber, the time-of-fiight system of scintil- 
lation counters, the particle identification system based 
on the aerogel Cherenkov counters, the calorimeter (the 
liquid krypton in the barrel part and the Csl crystals in 
the end caps) and the muon tube system inside and out- 
side of the magnet yoke. In this experiment the magnetic 
field was off and the liquid krypton calorimeter was out 
of operation. 

To suppress the machine background to the acceptable 
level, the following trigger conditions were used by OR 

1. signals from > 2 barrel scintillation counters coin- 
ciding with the Csl calorimeter signal, 

2. coinciding signals of two Csl end-caps, 

with the Csl energy threshold of about 75 MeV. The 
Monte Carlo simulation employing the JETSET-7.4 code 
|25] yields the trigger efficiency of about 0.4 for J/ip de- 
cays and about 0.43 for ip' decays. 

B. Multihadronic event selection 

For the off-line event selection the following conditions 
were applied: 

1. > 3 charged tracks or 2 acolinear charged 
tracks (cos 8 < 0.95) from the interaction region 
(p < 5 mm, \z\ < 120 mm), 

2. the Csl energy > 1.15 of the hardware threshold. 

The second condition serves to exclude the hardware 
threshold instability. The detection efficiency, deter- 
mined by the visible peak height and the table value 
of the leptonic width, is about 0.25 for J/ip (~ 20 • 10 3 
events) and about 0.28 for ip' (~ 6 • 10 3 events). 

The residual machine background (beam-gas and 
beam- wall) does not exceed 5 nb. The systematic er- 
ror in the meson masses related to its variation is less 



than 1 keV. To obtain this estimate, the background was 
increased by a few times by adding the appropriate frac- 
tion of unselected events to the selected ones at each 
experimental point (see Fig. |2J). Further suppression of 
the background leads to the detection efficiency loss and 
does not improve the mass accuracy. 

The meson mass value is rather sensitive to the de- 
tection efficiency variation during the energy scan, its 
reduction by 1% at one point causes the ip' mass shift 
up to 5 keV. To ensure the detection efficiency stability, 
all electronic channels having problems at any moment 
of the experiment were excluded from the off-line anal- 
ysis. Besides, the relative hit efficiencies of all detector 
subsystems were obtained for all experimental points us- 
ing the cosmic ray runs, the multihadron event statis- 
tics and (when possible) background events. These ef- 
ficiencies were applied to the real multihadron events to 
determine the relative point-specific correction factors. 
The variation of the drift chamber spatial resolution was 
handled similarly. 

The correction procedure described above shifts the 
mass by (+6.3 ± 2.3) keV for J/tp and (+0.2 ± 2.0) keV 
for ip' . The errors include statistical errors of the hit 
efficiency determination and the uncertainties of the cor- 
rection procedure employed. The shift of the J/ip mass is 
mainly due to the false alarm of the safety system which 
stopped the gas flow in the drift chamber (one point of 
the second scan). The values of shifts are given just for 
information, only the errors are of importance. 



C. Luminosity measurements 

For the operative VEPP-4M luminosity measurements 
single bremsstrahlung monitors were installed in both 
e + - and e~-directions. Their stability is not sufficient 
for the precision mass measurements, therefore Bhabha- 
scattering events detected by the end-cap Csl calorime- 
ter were employed. The fiducial polar angle interval is 
17.5° < 6 < 35°. The resonance contribution is not neg- 
ligible in this angular range, so that the e + e~ — ► ip — > 
e + e~ and e + e~ — * e + e~ interference should be taken into 
account. The correction to the number of e + e~-events 
has been included according to 26] . The values of the 
total and electronic widths were taken from [27j and the 
beam energy spreads were known from the experiment. 

Unlike e + e~ — * the e + e~ — > e + e~ interference 

dip is at the high-energy side of the resonance curve, 
so the mass shifts due to the correction are positive: 
(15 ± 1) keV (J/ip ) and (5 ± 0.5) keV (ip'). 

To the errors quoted, 2 keV should be added for J /ip 
and 3 keV for ip' to cover the calorimeter instabilities. 
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VII. RESONANCE EXCITATION CURVE 
FITTING 

A. Introduction 

Taking into account the beam energy spread, the event 
production rate for the collider energy E can be written 
as 



F(E) = Ja(W) dL(E,W) 



(10) 



where W is the cm. energy of the collision, dL(E, W) is 
the differential luminosity and a(W) is the cross section. 

In case of the narrow vector meson production in the 
reaction e + e~ — > V — * hadrons 



a{w) = wj dxx 



(W(l- x)- M) 2 +r 2 /4 



(ii) 



where T, r ee and are total and partial widths of the 
meson, M is its mass and &(x, W) is the probability to 
lose the fraction of energy x because of the initial state 
radiation [2^ (we substituted the Breit-Wigner cross sec- 
tion with the physical value of T ee including the vacuum 
polarization effects and used W instead of s = W 2 ). 

After the corrections introduced in Sec. E] to exclude 
the asymmetry causing the resonance mass shift, the 
symmetric expression for the differential luminosity can 
be used: 



dL(E,W) L(l + k(W-2Ef 



dW 



2%a w (I,J) 

(W-2E) 2 



exp 



2a 2 w (I, J) 



(12) 



where the (free) small parameter k is introduced to cover 
non-gaussian effects due to the /3-function chromaticity 
and other possible reasons. The collision energy spread 
aw can depend on the beam current I because of the mi- 
crowave instability reviewed in |29| and/or on the current 
density J due to the multiple intra-beam scattering (|3(j 
and references therein). The latter depends not only on 
the beam current, but also on the beam sizes modify- 
ing substantially by the collision effects, thus it must be 
considered as an independent parameter. 

Formula (|11|) ignores the interference between the reso- 
nant and nonresonant hadron production. With the suf- 
ficient accuracy [2(j: 



a(W) 



12tt 



(0 



TM 



Imf- 



2a\/RT ee T h 
Ml 



A -Re 



(13) 



Here a is the fine structure constant, R = a^'/a 



A denotes the fraction of events interfering with the non- 
resonant hadronic cross section and 



7r \ m P 2 



f 



M/2 



-W + M- iT/2 



1-/3 



In the limit of zero resonance width r — > 0, assuming 
the Gaussian energy spread and ignoring the interference, 
it is possible to express (|10|) in terms of known functions: 



F(E) 



6tt 2 r ee r h (tow\ p 
m 2 r \ m ) 

T{\ + [3) f (W-2E) 2 



2tt aw 



exp 



4^ 



W-2E 



aw 



(1 



where V is the gamma-function, D^p 
parabolic cylinder function and 



(14) 

-S)L, 

is the Weber 



S--1 — -1 

~ 7T I 3 ~ 2 



B. Interference effect treatment 

In principle, the interference magnitude can be left free 
in the fit and extracted from the data together with the 
resonance mass and the machine energy spread aw ■ Un- 
fortunately, this affects too much the statistical accuracy 
of the mass measurements, so one has to fix it. The zero 
magnitude is usually assumed. 

In this analysis we fix the interference parameter A, em- 
ploying the parton model of the onium decays. It assumes 
that J/ip decays to the light gg-pairs with the probability 
of RB^^ and to the gluon triplet ggg or the ggj mixture 



2)B^, where is the 



with the probability of 1 — [R 
muon branching ratio. The events J/?/ 1 — > qq are iden- 
tical to those in the non-resonant continuum, and there 
is the 100% interference in this case. For the hypotheti- 
cal heavy onium decays, ggg-qq interference is negligible 
due to the difference in the angular distributions (three 
jets vs two jets). For the real J/i(>, the angular distri- 
butions do not differ much but the interference phases 
are individual for all exclusive final states, therefore the 
net interference effect should be small just because of the 
large number of decay modes. 

So, the fraction of multihadron J/tp decay events inter- 
fering with the continuum A w RB^j (1 — 2B MM ) as 0.17. 
The uncertainty of A related to the final number of 
the decay modes was estimated by the multiple assign- 
ment of the arbitrary interference phases to all J/ip 
and decay modes implemented in the JETSET-7.4 
Monte Carlo code [25|. The following values have been 



obtained: A 



0.17 ±0.03 and A 



0.023 ±0.009. 



The corresponding mass shifts are (±7.0 ± 1.3) keV and 
(±2.0 ± 0.8) keV, respectively. The accuracy of the par- 
ton model predictions used gives a small addition to the 



9 



1 














1 


Scan IV: 

o w = om 


i 






1 




7 VT 


Scans 1,11: 
<i> = -0.011 
<j> = -0.186 




1 








°w = 

± 


0.839 
0.013 












'"' Scan 
i <i> 
V <j> 


III: 

= +0.013 












\\ : 


= +0.246 
= 0.900 
: 0.020 










, , , 


, , , 







1546 



1547 



1548 



1549 



1550 



1551 
E, MeV 



250 



200 



150 



100 



50 





























. Sea 


ns I-m: 














/ °w 


= 1.330 
± 0.024 






















...... f" 








................... 


................... 





1840 



1841 



1842 



1843 



1844 



1845 1846 
E, MeV 



FIG. 7: The results of four J/ip scans (the energy spread 
values ow and the mean i, j values for the first two and the 
third scans are presented). 



FIG. 8: The results of three ip' scans (the mean energy spread 
value ow is presented). 



errors quoted. The fitting procedure automatically shifts 
the mass value so only the errors of the quoted values are 
of importance. 



C. On collision energy spread variations 

The single beam measurements of the longitudinal and 
radial bunch sizes indicate that the energy spread in 
VEPP-4M depends on the bunch current. In the J/ip re- 
gion the dimensionless slope (da e / 'a e) / '(dl / (I)) — 0.07. 
As the current decreases during the run, it leads to the 
symmetric distortion of the collision energy distribution 
and contributes to the k parameter of (|12fl . If the mean 
value of the energy spread is not the same at all energy 
points where data were collected, the fake mass shift can 
appear. 

To take into account the energy spread variations in 
the resonance curve fit, we assumed that it linearly de- 
pends on the beam current and the current density in the 
vicinity of their mean values 



a E k. (a E ) (l + Oi-i + Oij-j), 

1 i J 

1 = 1—, 1. 7 = ■; — ; 

(I) ' (J) 



(15) 



The product I ■ £ n with the specific luminosity 
£ = L/I+I- oc l/cr* cr* were chosen as a measure of the 
current density effect as the beam sizes in the ring were 
not permanently monitored. The cti and ay were con- 
sidered as free parameters, the values n = 2 and n — 1 
were tried. If the synchrotron contribution in the radial 
size dominates, and the vertical beam size is due to the 
coupling (cr* oc a*p), one has n = 2 for the multiple 
intra-beam scattering. With n = 1 the j parameter just 
characterizes the strength of the collision effects (I-l oc £ y 
for the flat beams). 



D. Fitting procedure 

Each data acquisition run was subdivided into subruns 
with a minor variation of the beam currents. For each 
subrun the beam energy E was assigned and the parame- 
ters i± and j± were calculated for the e + - and e~-beams. 

The observed number of the multihadron events N 
has been fitted as the function of E, i and/or j using 
the log-likelihood method. The calculations of the ex- 
pected number of the resonance events were performed 
by the numerical convolution of l|13fl an d (|12fl with 



&W — y 17 e+ + °E- depending on i± and/or j± accord- 
ing to l|15fl . The free parameters were: the constant con- 
tinuum cross section a c , the event detection efficiency e 
for the given value of the leptonic width T ee , the reso- 
nance mass M and the energy spread parameters (<7w), 
on and/or a>j. The PDG table values |27j were used for 
the resonance widths T and T ee . The A parameter in 113|) 
was fixed as described above. 

For the fitting procedure verification and the system- 
atic error checks, N(E) fits were performed using Ijl4|l 
with the results presented in Fig. [7| and Fig. [SJ 

The peak height for the third J/tp scan differs from 
those for the first and the second ones due to the aw 
dependence on the parameter j. The energy spread in 
the fourth scan was decreased intentionally (see Sec. IIIip . 

In the tp' case, the energy spread variations are not 
seen due to higher energy and narrower i- 1 j-ranges. 

The chi-squared values of the fits are satisfactory 
(P(x 2 ) > 0-1) f° r a U J/i^ an( i i 3 ' scans even if the depen- 
dence of the energy spread on the beam current within a 
scan is ignored. 

The detection efficiencies obtained by the fit assum- 
ing the world average values of the leptonic widths 



= 5.26 ±0.37 keV, T 



2.19 ±0.15 keV W 



agree with those obtained by Monte Carlo simulation 
within their errors. The systematic error of the Monte 
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FIG. 9: J/ip mass values for each scan relative to the world 
average value |2?j ]. the resulting average values are shown by 
the dashed lines (statistical errors only). 



FIG. 10: ip' mass values for each scan relative to the world 
average value [2^ |. the resulting average values are shown by 
the dashed lines (statistical errors only). 



Carlo calculations and the error of the absolute luminos- 
ity calibrations in this experiment are large (about 12% 
in total), so we do not present our leptonic width values. 
For the mass measurements these errors are not impor- 
tant. 



VIII. MEASURED MASS VALUES AND ERROR 
DISCUSSION 

The mass values obtained in each scan assuming the 
constant value of aw during a scan and the resulting 
average mass values for J/ip and ip' as well as their sta- 
tistical errors are presented in Figs. 191101 

To obtain the resulting averaged mass values, the scans 
were considered as independent experiments. The indi- 
vidual mass values of the scans were weighted using their 
statistical errors and ignoring the systematic ones. Corre- 
spondingly the systematic errors of the individual scans 
were weighted. Such procedure overestimates the total 
error, but allows one to separate the statistical and sys- 
tematic errors of the resulting value. 

The mass values for all scans are in good agreement 
even when the difference in systematic errors is ignored 
(the x 2 is given in Figs. l9ll'D|l . The resulting statistical 
accuracy is 10 keV for J/ip and 25 keV for ip' . 

The fit of J/ip scans I-^III does not show 
the statistically significant direct dependence on / 
(a, = 0.037 ± 0.055 ) and gives the dimcnsionless slope 
a, =0.059 ± 0.018 similar to that of the single-beam 
measurements. The mass deviation from the value for 
ai,aj = does not exceed 1.5 keV. The correction is 
rather small, so the mass values for a^, ay = have been 
assigned to the scans I-III of J/ip with the systematic 
error of 1.5 keV related to the energy spread variations. 



For the J '/tp scan IV with the damping decrements re- 
arranged, one has ctj = 0.45 ± 0.15 for at fixed at zero and 
cti = —0.13 ± 0.16 for otj fixed at zero with the mass vari- 
ations of +10 keV and —8 keV, respectively. The mass 
value for ai ,aj — has been assigned to the scan IV of 
J/ip with the systematic error of 10 keV related to the 
energy spread variations. The mass value obtained as- 
suming the dependence on J oc I- 1 1 does not contradict 
to this error estimate. The difference of the J/ip mass 
values obtained in the scan IV and in the scans I-^III is 
(-11 ±22 ±12) keV. 

The systematic uncertainties are listed in Table IIIII 
The sources of uncertainties not exceeding 0.3 keV (the 
uncertainties of the world average values of particle prop- 
erties used, the resonance mass uncertainty related to (@J, 
the radiative corrections uncertainties etc.) are omitted. 

The weighted contribution of the energy spread vari- 
ation to the systematic error of the resulting J/ip mass 
value is about 3 keV. During the ip' scans, the statisti- 
cally significant variations of the energy spread have not 
been observed, a relatively big systematic error of 2 keV 
is due to the bigger value of aw ■ 

The uncertainty of the energy assignment includes the 
statistical component which was estimated using the mul- 
tiple data fits with a randomly generated energy devia- 
tion. The systematic contribution consists of two parts. 
The first one is related to the energy prediction function 
(p.f.) choice and the energy assignment policy. Three 
values of energy were tried for each run: the unbiased p.f. 
value and two values of the shifted p.f. which exactly re- 
produce the preceding and following calibration results. 
Besides, the best p.f. parameter variation within their 
errors were allowed and a few different functions were 
tried. The second part is due to the difference between 
the data acquisition runs and the energy calibration runs 
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TABLE III: The systematic uncertainties in J/ip and ip' masses (keV) 



Source 


J/xl> 


i>' 


Energy spread variation f Sec. IVIIII IVII CI 


3.0 


2.0 


Encrev assignment: statistical uncertainty (Sec. IVIIII IIV Bll 


2.5 


3.5 


Enerev assignment: prediction function choice (Sec. 1 Villi IIV Bl 


2.7 


1 -J 
1. 1 


Energy assignment: radial betatron oscillations (Sec. IIV Tit 


<1.5 


<1.8 


Energy assignment: beam separation in the additional LP. (Sec. ILVFH 0.4* 


0.4* 


Beam misalignment in the interaction point (Sec. fvl 


1.8 


5.1 


e~*~-, 6 -energy difference (Sec. IVHH 


o n 
< z.U 


<2.0 


Non-gaussian collision energy distribution (Sec. IVIIII IVII Afl 


<1.5 


<2.0 


/9-function chromaticity (Sec. IvTl 


2.0* 


2.5* 


Beam potential (Sec.lVt 


1.0* 


1.0* 


Single energy calibration (Sec. Ill Cll 


0.6 


0.8 


Detection efficiency instability (Sec. IVIFil 


2.3 


Z.U 


Luminosity measurements (Sec.lVl'Ul 


2.2 


3.0 


Interference in the hadronic channel (Sec. IVII a\ 


1.3 


0.8 


Hesidual machine background (Sec. IVI Bt 


<1.0 


< 1.0 


Sum in quadrature 


f«7.3 


« 8.9 


* — correction uncertainty 






TABLE IV: The corrections applied to the J/ip and ip' mass 


values (keV). 




Correction for J/ip(l~lIl) 


j/V(iv) 


if 


Vertical orbit disturbances (Sec. Ill CI — 0.8 ± 0.6 - 


-0.8 ±0.6 -0.6 ±0.4 


Separation in the additional LP. (Sec. IIV Bt -3.4 ± 0.4 - 


-3.4 ±0.4 -4.0 ±0.4 


3- function chromaticity (Sec.fvTl -4.0 ± 2.0 - 


-1.5 ±0.7 5.0 ±2.5 


Beam potential (Sec.EJ 1.9 ± 1.0 


2.1 ±1.0 2.0 ±1.0 


Total -6.3 ±2.4 - 


-3.6 ±1.7 1.4 ±2.8 



(the radial betatron oscillations and the separation in the 
additional LP., Sec. llVB)) . 

The mass uncertainty caused by the beam misalign- 
ment in the experimental LP. and the electrostatic dis- 
persion was evaluated similarly to that of the statistical 
error of the energy assignment. 

Estimates show that the difference of the energies of 
e + - and e~ -beams is small. The value shown in Ta- 
ble [TH] has been obtained using a few energy calibrations 
with the e + -bcam performed during the stability runs 
(Sec. llVAl 

The symmetric distortion of the distribution in the col- 
lision energy W can shift the mass due to the asymmetry 
of radiative corrections. The corresponding uncertainty 
has been estimated by leaving the k parameter of H12|) 
free (the values shown in Fig. 191101 are for k = 0). 

All other sources of the systematic error were discussed 
above. The resulting systematic error in the mass is 
about 7 keV for J/tp and about 9 keV for tp'. The correc- 
tions applied to the fit results are presented in Table llVl 

The results obtained can be presented in the form 

Mj^ - M^ G = 47 ± 10 ± 7 (±40 ) keV, 
M r - M^P G = 151 ± 25 ± 9 (±90 ) keV, 

demonstrating the agreement with the world average val- 
ues. 



IX. CONCLUSION 

The new high precision measurement of the J/ip- 
and '(//-meson masses has been performed at the collider 
VEPP-4M using the KEDR detector. The following mass 
values have been obtained: 

M J/4 , = 3096.917 ± 0.010 ± 0.007 MeV 
Afy = 3686.111 ±0.025 ±0.009 MeV 

The relative measurement accuracy reached 4 • 10~ 6 for 
J/ip , 7 ■ 10~ 6 for if}' and is approximately 3 times better 
than that of the previous precise experiments. 
For the mass difference our result is 

Afy - M J/4 , = 589.194 ± 0.027 ± 0.011 MeV. 
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